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Abstract 



^ \ We carry out a high-precision simulation of the two-dimensional SU(3) prin- 

cipal chiral model at correlation lengths £ up to ~4 x 10 5 , using a multi-grid 
Monte Carlo (MGMC) algorithm. We extrapolate the finite-volume Monte 
Carlo data to infinite volume using finite-size-scaling theory, and we discuss 
carefully the systematic and statistical errors in this extrapolation. We then 
compare the extrapolated data to the renormalization-group predictions. For 
£ > 10 3 we observe good asymptotic scaling in the bare coupling; at £ ~ 4 x 10 5 
the nonperturbative constant is within 2-3% of its predicted limiting value. 



PACS number(s): ll.10.Gh, 11.15.Ha, 12.38.Gc, 05.70.Jk 



A key tenet of modern elementary-particle physics is the asymptotic freedom 
of four-dimensional nonabelian gauge theories 0. However, the nonperturbative 
validity of asymptotic freedom has been questioned @; and numerical studies of lattice 
gauge theory have thus far failed to detect asymptotic scaling in the bare coupling 
||. Even in the simpler case of two-dimensional nonlinear cx-models [|J, numerical 
simulations at correlation lengths £ ~ 10-100 have often shown discrepancies of order 
10-50% from asymptotic scaling. In a recent paper || we employed a finite-size- 
scaling extrapolation method || [7|, [|, || to carry simulations in the 0(3) a-model to 
correlation lengths £ ~ 10 5 ; the discrepancy from asymptotic scaling decreased from 
~ 25% to ~ 4%. In the present Letter we apply a similar technique to the SU(3) 
principal chiral model, reaching correlation lengths £ ~ 4 x 10 5 with errors < 2%. 
For £ > 10 3 we observe good asymptotic scaling in the bare parameter /3; moreover, 
at £ ~ 4 x 10 5 the nonperturbative ratio ^observed/ ' £,theor,3-ioo P is within 2-3% of the 
predicted limiting value. 

We study the lattice cx-model taking values in the group SU(N), with nearest- 
neighbor action 7~C{U) = — /3X}Retr([/J U y ). Perturbative renormalization-group 
computations predict that the infinite- volume correlation lengths £( exp ) and £^ [ 10 1 
behave as 



v N J 

as P — > oo. Three- loop perturbation theory yields 

ai = -0.121019iV + 0.725848iV" 1 - 1.178097iV _3 . (2) 
The nonperturbative constant C^ V ) has been computed using the thermodynamic 



Bethe Ansatz [13]: 



y/e n/N f 
16v/tF sin(7r/iV) 6XP \ ~2iV 2 



C &*rt = Tc~~/z: „:„/_/An ex P( ~ n oa^ I ' ( 3 ) 



The nonperturbative constant C*(2) is unknown, but Monte Carlo studies indicate 
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that C^/C^ lies between « 0.985 and 1 for all iV > 2 Q; for iV = 3 it is 
0.987 ± 0.002 pf. Monte Carlo studies |J, 0, of the SC/(3) model up to 



^ ~ 35 have failed to observe asymptotic scaling ([1]); the discrepancy from (|1])-([|) is 
of order 10-20%. 

Our extrapolation method is based on the finite-size-scaling (FSS) Ansatz 

= F (t{^L)/L- s) + 0(C",L-) , (4) 

where O is any long-distance observable, s is a fixed scale factor (here s = 2), L 
is the linear lattice size, Fq is a universal function, and u is a correction-to-scaling 
exponent. We make Monte Carlo runs at numerous pairs (/3, L) and (/3, sL); we then 
plot O(0, sL)/0(/3, L) versus L)/L, using those points satisfying both L) > 
some value £ m j n and L > some value L min . If all these points fall with good accuracy 
on a single curve, we choose a smooth fitting function Fq. Then, using the functions 
F% and F , we extrapolate the pair (£, C) successively from L — > sL — > s 2 L 
oo. See [§] for how to calculate statistical error bars on the extrapolated values. 
We have chosen to use functions Fq of the form 

F {x) = 1 + a x e- x l x + a 2 e- 2 ' x + ... + a n e~ n l x . (5) 

We increase n until the x 2 of the fit becomes essentially constant; the resulting x 2 
value provides a check on the systematic errors arising from corrections to scaling 
and/or from inadequacies of the form ([|). The discrepancies between the extrapolated 
values from different L at the same f3 can also be subjected to a x 2 test. Further 
details on the method can be found in || || . 

We simulated the two-dimensional SU(3) cx-model using an XF-embedding multi- 
grid Monte Carlo (MGMC) algorithm @. We ran on lattices L = 8,16,32,64, 



128, 256 at 184 different pairs (/3, L) in the range 1.65 < j3 < 4.35 (corresponding to 
5 ^ ^oo ~ 4 x 10 5 ). Each run was between 4 x 10 5 and 5 x 10 6 iterations, and the 



total CPU time was one year on a Cray C-90 p0[. The raw data will appear in [21 



Our FSS data cover the range 0.08 < x = / L < 1.12, and we found ten- 
tatively that for O = £ a thirteenth-order fit (|5|) is indicated: see Table |I[ There 
are significant corrections to scaling in the regions x ^ 0.84 (resp. 0.64, 0.52, 0.14) 
when L = 8 (resp. 16, 32, 64): see the deviations plotted in Figure |l[ We therefore 
investigated systematically the x 2 of the fits, allowing different cuts in x for differ- 
ent values of L: see again Table |I|. A reasonable x 2 is obtained when n > 13 and 
Xmin > (0.80, 0.70, 0.60, 0.14, 0) for L = (8, 16, 32, 64, 128). Our preferred fit is n = 13 
and x min = (oo, 0.90, 0.65, 0.14, 0): see Figure @, where we compare also with the 
perturbative prediction 



(6) 



valid for x > 1, where a = 2N/(N 2 - 1), w = N/(8n) and w x = N 2 /(128tt 2 ). 

The extrapolated values £jg from different lattice sizes at the same (3 are consistent 
within statistical errors: only one of the 58 (3 values has a x 2 too large at the 5% 
level; and summing all f3 values we have x 2 = 64.28 (103 DF, level = 99.9%). 

In Table § we show the extrapolated values £w from our preferred fit and some 
alternative fits. The deviations between the different fits (if larger than the statisti- 
cal errors) can serve as a rough estimate of the remaining systematic errors due to 
corrections to scaling. The statistical errors in our preferred fit are of order 0.5% 
(resp. 0.9%, 1.1%, 1.3%, 1.5%) at w 10 2 (resp. 10 3 , 10 4 , 10 5 , 4 x 10 5 ), and the 
systematic errors are of the same order or smaller. The statistical errors at different 
j3 are strongly positively correlated. 

In Figure | (points + and x) we plot CJ, est imate (00,0.90,0.65) divided by the two- 
loop and three- loop predictions (0)~(0) for ^ exp \ The discrepancy from three- loop 
asymptotic scaling, which is ~ 13% at (3 = 2.0 (^ ~ 25), decreases to 2-3% at (3 = 
4.35 (£oo ~ 3.7x 10 5 ). For (3 ^ 2.2 (^ > 60) our data are consistent with convergence 
to a limiting value C^/C^xp) ~ 0.99-1 with the expected 1//3 2 corrections. 
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We can also try an "improved expansion parameter" [^2, [l2j based on the energy 
E = iV _1 (Retr(£/o First we invert the perturbative expansion |T3 



E{(3) = 1- 



N 2 -l 

4N/3 



N 2 - 2 0.0756 — 0.0634 A^ 2 + 0.01743 iV 4 , 

1 + w + wjr + 0{1/p > 



(7) 

and substitute into ([I]); this gives a prediction for £ as a function of 1 — £7. For _E 
we use the value measured on the largest lattice (which is usually L = 128); the 
statistical errors and finite-size corrections on E are less than 5 x 10~ 4 , and they 
induce an error less than 0.85% on the predicted ^ (less than 0.55% for (3 > 2.2). 
The corresponding observed/predicted ratios are also shown in Figure ^ (points □ 
and O). The "improved" 3-loop prediction is extremely flat, and again indicates a 
limiting value ~ 0.99. 

Further discussion of the conceptual basis of our analysis can be found in pi. 
Details of this work, including an analysis of the susceptibility \i will appear elsewhere 

& 
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Figure 1: Deviation of points from fit to with s — 2, n — 13, x m j„ = 
(oo,oo,oo,0.14,0). Symbols indicate L = 8 (F^), 16 (*#x), 32 (+). Error bars 
are one standard deviation. Curves near zero indicate statistical error bars (± one 
standard deviation) on the function F^(x). 
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Figure 2: £(/3, 2L) /£(/?, L) versus £(J3,L)/L. Symbols indicate L = 8 ffl), 16 
(%£x), 32 (+), 64 (x), 128 (□). Error bars are one standard deviation. Solid 
curve is a thirteenth-order fit in (|5|), with x m i n = (oo, 0.90, 0.65, 0.14, 0) for L = 
(8, 16, 32, 64, 128). Dashed curve is the perturbative prediction (^). 
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Figure 3: ^ Mtimate(o0) o.oo 1 o.6B)/^£S»r versus P- Error bars are one standard devia- 
tion (statistical error only). There are four versions of ^^ eor : standard perturbation 
theory in 1/(3 gives points + (2-loop) and x (3- loop); "improved" perturbation the- 
ory in 1 — E gives points □ (2-loop) and O (3-loop). Dotted line is the Monte Carlo 
prediction C^)/C^ p) = 0.987 ±0.002 [12[. 
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%min 


n = 11 


n = 12 


n = 13 


n = 14 


n = 15 


(0.50,0.40,0) 


180 718.80 
3.99 0.0% 


179 626.60 
3.50 0.0% 


178 560.20 
3.15 0.0% 


177 558.60 
3.16 0.0% 


176 558.30 
3.17 0.0% 


(oo,0.40,0) 


154 673.80 
4.38 0.0% 


153 566.30 
3.70 0.0% 


152 533.00 
3.51 0.0% 


151 532.10 
3.52 0.0% 


150 531.80 
3.55 0.0% 


(oo,oo,0) 


108 236.00 
2.19 0.0% 


107 172.40 
1.61 0.0% 


106 154.80 
1.46 0.1% 


105 154.70 
1.47 0.1% 


104 153.40 
1.48 0.1% 


(0.70,0.55,0.45) 


162 288.30 
1.78 0.0% 


161 219.20 
1.36 0.2% 


160 183.00 
1.14 10.3% 


159 182.50 
1.15 9.8% 


158 182.30 
1.15 9.0% 


(0.75,0.60,0.50) 


150 222.40 
1.48 0.0% 


149 172.20 
1.16 9.4% 


148 129.90 
0.88 85.6% 


147 129.80 
0.88 84.3% 


146 129.80 
0.89 82.9% 


(0.80,0.70,0.60) 


129 173.90 
1.35 0.5% 


128 135.00 
1.05 32.0% 


127 96.30 
0.76 98.1% 


126 96.28 
0.76 97.7% 


125 94.31 
0.75 98.1% 


(0.95,0.85,0.60) 


111 150.30 
1.35 0.8% 


110 107.20 
0.97 55.8% 


109 77.62 
0.71 99.0% 


108 77.62 
0.72 98.8% 


107 75.67 
0.71 99.1% 


(1.00,0.90,0.60) 


105 139.20 
1.33 1.4% 


104 100.90 
0.97 56.7% 


103 70.74 
0.69 99.4% 


102 70.73 
0.69 99.2% 


101 67.50 
0.67 99.6% 


(oo,0.90,0.65) 


92 130.00 
1.41 0.6% 


91 77.01 
0.85 85.2% 


90 60.85 
0.68 99.2% 


89 58.66 
0.66 99.5% 


88 58.31 
0.66 99.4% 


(00, 00, 0.65) 


78 96.09 
1.23 8.1% 


77 56.51 
0.73 96.2% 


76 49.55 
0.65 99.2% 


75 46.63 
0.62 99.6% 


74 45.94 
0.62 99.6% 


(00,00,00) 


52 55.85 
1.07 33.2% 


51 25.23 
0.49 99.9% 


50 25.17 
0.50 99.9% 


49 24.11 
0.49 99.9% 


48 24.10 
0.50 99.8% 



Table 1: Degrees of freedom (DF), x 2 , % 2 /DF and confidence level for the n th - 
order fit (|]) of £(/?, 2L) /£(/?, L) versus £({3,L)/L. The indicated x m i n values apply 
to L — 8,16,32, respectively; we always take x m i n = 0.14,0 for L = 64,128. Our 
preferred fit is shown in italics; other good fits are shown in sans-serif; bad fits are 
shown in roman. 
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771 in 


1.80 


2.00 


2.20 


2.40 


2.60 


2.85 


3.00 


3.15 


(0.70,0.55,0.45) 


10.455 (0.022) 


24.903 (0.066) 


57.13 (0.17) 


129.68 (0.41) 


290.5 (1.0) 


794.9 (3.2) 


1460 ( 6) 


2687 (11) 


(0.75,0.60,0.50) 


10.454 (0.022) 


24.886 (0.071) 


57.50 (0.18) 


130.83 (0.43) 


293.0 (1.1) 


801.7 (3.4) 


1473 ( 6) 


2709 (12) 


(0.80,0.70,0.60) 


10.450 (0.021) 


24.875 (0.073) 


57.41 (0.22) 


130.93 (0.64) 


293.6 (1.6) 


805.9 (5.0) 


1482 ( 9) 


2727 (17) 


(0.95,0.85,0.60) 


10.451 (0.021) 


24.870 (0.071) 


57.40 (0.21) 


130.93 (0.63) 


293.7 (1.6) 


806.6 (6.1) 


1483 (12) 


2749 (25) 


(1.00,0.90,0.60) 


10.450 (0.022) 


24.872 (0.069) 


57.40 (0.21) 


130.94 (0.63) 


293.6 (1.6) 


806.8 (5.9) 


1484 (12) 


2749 (25) 


(oo,0.90,0.65) 


10.446 (0.022) 


24.859 (0.072) 


57.40 (0.21) 


131.00 (0.66) 


295.2 (2.1) 


809.6 (6.7) 


1489 (13) 


2761 (27) 


(oo, oo,0. 65) 


10.447 (0.022) 


24.863 (0.074) 


57.40 (0.22) 


131.01 (0.66) 


295.0 (2.1) 


809.7 (6.9) 


1487 (14) 


2759 (28) 


(00,00,00) 


10.454 (0.022) 


24.881 (0.074) 


57.39 (0.22) 


130.78 (0.66) 


295.6 (2.3) 


812.7 (9.8) 


1482 (22) 


2777 (49) 





3.30 


3.45 


3.60 


3.75 


3.90 


4.05 


4.20 


4.35 


(0.70,0.55,0.45) 


4957 ( 23) 


9117 ( 46) 


16780 ( 92) 


30959 (182) 


56766 ( 362) 


105205 ( 707) 


196197 (1396) 


360864 ( 2792) 


(0.75,0.60,0.50) 


4995 ( 24) 


9199 ( 47) 


16938 ( 93) 


31258 (185) 


57265 ( 366) 


106093 ( 736) 


197949 (1419) 


363905 ( 2880) 


(0.80,0.70,0.60) 


5032 ( 32) 


9268 ( 62) 


17066 (118) 


31492 (239) 


57687 ( 456) 


106807 ( 878) 


199117 (1690) 


366159 ( 3309) 


(0.95,0.85,0.60) 


5109 ( 49) 


9411 ( 92) 


17359 (178) 


32059 (346) 


58748 ( 650) 


108781 (1237) 


202868 (2360) 


372553 ( 4392) 


(1.0,0.90,0.60) 


5110 ( 51) 


9365 ( 99) 


17299 (196) 


31816 (372) 


58308 ( 702) 


107789 (1312) 


200994 (2493) 


369579 ( 4697) 


(oo,0.90,0.65) 


5132 ( 55) 


94 07 (105) 


17377 (208) 


31908 ( 398) 


58594 ( 766) 


108952 (1452) 


201796 ( 2817) 


371706 ( 5457) 


(oo,oo,0.65) 


5125 ( 55) 


9391 (110) 


17389 (229) 


32008 (463) 


58804 ( 941) 


109440 (1886) 


204587 (3779) 


376704 ( 7722) 


(00,00,00) 


5063 (102) 


9295 (217) 


16991 (447) 


30912 (903) 


55976 (1828) 


104740 (3678) 


192664 (7358) 


359299 (14787) 



Table 2: Estimated correlation lengths ^ as a function of (3, from various extrapola- 
tions. Error bar is one standard deviation (statistical errors only). All extrapolations 
use s = 2 and n = 13. The indicated x min values apply to L — 8, 16, 32, respectively; 
we always take x m i n = 0.14,0 for L = 64, 128. Our preferred fit is shown in italic; 
other good fits are shown in sans-serif; bad fits are shown in roman. 
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